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Abstract — The  capacity  optimization  problem  of  MIMO  links 
with  interference  has  attracted  an  increasing  interest.  Due  to 
the  nonconvexity  of  the  capacity  problem,  only  suboptimal 
solutions  can  be  found.  In  the  previous  works,  a  Gradient 
Projection  (GP)  algorithm  [1]  and  a  Quasi-Newton  (QN) 
method  [2]  were  proposed  to  provide  suboptimal  solutions 
subject  to  the  constant  power  constraint. 

In  this  paper,  we  derive  the  capacity  for  MIMO  links 
suffered  from  cochannel  interference  where  each  MIMO  link  is 
decomposed  via  SVD.  Then,  each  eigenchannel  of  MIMO  link 
is  represented  by  a  set  of  logical  links  with  a  set  of  discrete 
data  rates  and  discrete  powers.  An  Integer  Programming  based 
algorithm  (named  as  IP)  is  presented  to  solve  the  capacity 
optimization  problem.  The  solution  specifies  the  set  of  logical 
links  that  can  transmit  simultaneously. 

Numerical  results  show  that  GP  and  QN  methods  achieve 
better  performance  than  IP  method  for  the  case  of  weak 
interference  because  of  the  convexity  of  the  optimization 
problem  when  INR  is  sufficiently  small.  In  the  case  of  strong 
interference,  IP  method  achieves  better  performance  than  GP 
and  QN  methods,  which  means  that  transmitting  one  link  at  a 
time  is  better  than  transmitting  all  links  simultaneously  with 
full  power.  In  other  words,  scheduling  links  to  transmit  is  more 
efficient  for  the  case  of  strong  interference. 

I.  Introduction 

Multiple-Input  Multiple-Output  (MIMO)  systems  have 
shown  great  promise  in  providing  high  spectral  efficiency 
for  single  user  wireless  link  without  interference  [3],  [4]. 
There  has  also  much  work  on  the  MIMO-based  cellular 
networks,  which  include  MIMO  multiple  access  (MIMO- 
MAC)  [5],  [6]  and  MIMO  broadcast  systems  [7],  [8].  Both 
systems  have  one  common  end  of  the  communication  link  - 
either  the  transmitter  of  MIMO-BS  or  the  receiver  of  MIMO- 
MAC.  There  has  been  great  interest  in  extending  the  MIMO 
communication  to  the  multi-user  systems  with  interference. 
The  transmission  scheme  of  each  user  depends  on  that  of 
other  users  since  the  interferences  at  each  user  depend  on 
all  the  transmit  covariance  matrices  in  the  network. 

We  consider  a  scenario  of  MIMO  system  where  the 
nodes  aim  to  maximize  the  total  system  capacity  rather  than 
the  individual  capacity.  It  is  well  known  that  total  system 
capacity  problem  appeared  in  the  MIMO-based  ad  hoc 
network  is  a  non-convex  programming  problem,  and  only 
suboptimal  solution  can  be  found.  A  Gradient  Projection 
(GP)  algorithm  [1]  and  a  Quasi-Newton  (QN)  method  [2] 
were  proposed  to  provide  suboptimal  solutions  subject  to 
the  constant  power  constraint.  However,  there  exist  some 
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special  cases  where  the  optimal  solutions  can  be  achieved. 
With  the  assumption  of  constant  power  constraint  at  each 
node.  Ye  and  Blum  [1]  show  that  the  total  system  capacity 
is  a  convex  function  of  links’  covariance  matrices  when 
the  interferences  from  other  links  are  sufficiently  small,  and 
hence  the  global  optimal  can  be  achieved.  They  also  show 
that  the  total  system  capacity  can  be  maximized  by  enforcing 
all  users  to  perform  beamforming  along  a  certain  direction 
when  the  interferences  from  other  links  are  sufficiently  large. 

In  this  paper,  we  derive  the  capacity  for  MIMO  links 
suffered  from  cochannel  interference  where  each  MIMO 
link  is  decomposed  via  SVD.  Then,  we  propose  an  Integer 
Programming  (IP)  method  to  study  the  total  system  capacity 
problem  for  the  MIMO  system  with  interference.  Both  trans¬ 
mitters  and  receivers  are  assumed  to  have  perfect  knowledge 
of  channel  state  information  (CSI).  Each  MIMO  link  is 
decomposed  into  parallel  independent  eigenchannel.  Then, 
each  eigenchannel  is  represented  by  a  set  of  logical  links, 
and  each  logical  link  is  determined  by  different  discrete 
powers  and  discrete  data  rates.  Also,  the  transmitters  of  each 
MIMO  link  are  subject  to  the  total  power  constraint.  The 
solution  specifies  the  set  of  logical  links  that  can  transmit 
simultaneously. 

The  remainder  of  this  work  is  organized  as  follows. 
Section  II  briefly  introduces  some  related  work.  In  Section 
III,  notation  and  problem  definition  are  given.  The  theoretical 
capacity  and  an  integer  programming  based  algorithm  to 
maximize  the  total  system  capacity  are  presented.  Then, 
Section  IV  shows  the  results  from  numerical  experiments. 
Finally,  concluding  remarks  are  given  in  Section  V. 

The  notations  in  this  paper  are  as  follows.  The  boldface 
denotes  matrices  and  vectors.  For  a  matrix  A,  Ar  denotes 
the  conjugate-transpose.  |A|  denotes  the  determinant,  and 
Tr( A)  denotes  the  trace  of  a  square  matrix.  AAO  means 
that  A  is  a  positive  semi-definite  matrix.  A  is  a  Hermitian 
matrix  if  A  =  AL  I  denotes  the  identity  matrix  with 
the  appropriate  dimension  from  the  context.  E[-]  denotes 
the  statistical  expectation.  C(x  x  y)  denotes  the  complex 
space  of  x  x  y  matrix.  min{x,y)  and  max(x,y)  denotes 
the  minimum  and  the  maximum  of  two  real  numbers  x  and 
y,  respectively. 

II.  Related  Work 

There  is  increasing  interest  in  the  MIMO-based  ad-hoc 
network.  In  [9],  Blum  investigates  the  MIMO  capacity  with 
interference  where  single-user  detection  is  assumed  at  the 
receiver.  Without  CSI  at  the  transmitter,  Blum  shows  that  the 
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optimum  is  either  the  optimum  interference-free  approach, 
which  puts  equal  power  into  each  antenna,  or  a  singular 
mode,  which  puts  all  powers  into  a  single  antenna,  if  the 
interference  is  either  sufficiently  weak  or  sufficiently  strong. 

In  [1],  Ye  and  Blum  study  the  asymptotic  behavior  of  the 
system  mutual  information  and  propose  a  Gradient  Projec¬ 
tion  (GP)  algorithm  to  maximize  the  total  system  capacity 
subject  to  the  constant  power  constraint.  In  [2],  a  Quasi- 
Newton  (QN)  method  is  proposed  to  solve  the  capacity 
problem  by  approximating  the  inverse  of  the  Hessian  matrix 
instead  of  computing  the  real  one.  QN  method  can  achieve 
the  super  linear  convergence  and  reduce  the  computation 
complexity.  However,  only  the  suboptimal  solutions  can  be 
found  due  to  the  non-convex  nature  of  the  optimization  prob¬ 
lem.  Furthermore,  the  constant  power  constraint  reduces  the 
achievable  maximum  system  capacity  in  the  case  of  strong 
interference.  In  [10],  Chen  and  Gains  study  the  asymp¬ 
totic  spectral  efficiency  of  K  simultaneously  communicating 
transmiter-receiver  pairs.  Without  CSI  at  the  transmitters, 
the  asymptotic  network  spectral  efficiency  is  limited  by  nr 
nats/s/Hz,  while  with  CSI  at  the  transmitters,  the  asymptotic 
network  spectral  efficiency  is  at  least  nt  +  nr  +  2fnt  +  nr 
nats/s/Hz. 

III.  System  Model  and  IP  Algorithm 
A.  MIMO  Link  with  Interferences 

In  this  paper,  we  consider  a  MIMO  interference  system 
with  L  links  where  each  node  is  equipped  with  nt  trans¬ 
mit  antennas  and  nr  receive  antennas.  Suppose  all  MIMO 
nodes  communicate  in  the  same  channel,  and  cochannel 
interference  is  presented  when  more  than  one  link  transmit 
simultaneously.  We  use  H.,  j  e  C(nr  x  nt)  to  represent  the 
channel  gain  matrix  from  the  transmit  antennas  of  link  j  to 
the  receive  antennas  of  link  i.  Assume  the  CSI  (or  HQQ 
corresponding  to  link  i  is  available  at  both  the  transmitter 
and  the  receiver  of  link  i. 

The  complex  base-band  signal  vector  received  by  the 
receiver  node  of  link  i  is  given  by 

L 

y i  =  \  “H,, x,  +  Y  +  ni i 

where  x,  represents  the  normalized  transmitted  signal  of  link 
i,  pi  denotes  the  signal-to-noise  ratio  (SNR)  of  link  i,  pi  j 
denotes  the  interference-to-noise  ratio  (INR)  of  link  i  due 
to  the  interference  from  link  j,  and  n,  is  the  noise  vector. 
The  entries  of  H,  j  and  n,  are  independent  and  identically 
distributed  complex  Gaussian  random  variables  with  zero 
mean  and  unit  variance.  For  simplicity,  we  assume  all  of  the 
interference  signals  x;  are  unknown  to  link  i. 

Let  Q,  be  the  covariance  matrix  for  the  transmit  signal 
vector  Xj,  and  defined  as 

Q;  =  -B[xjxf]. 

Tr(Qj)  specifies  the  transmit  powers  of  transmit  antennas 
of  link  i,  and  we  have  Q.,V0  and  Tr( QQ  <  1  where  1  rep¬ 


resents  the  normalized  maximum  power.  The  interference- 
plus-noise  of  link  i  is  x  /i,  ,IL  ,x,  +  n,,  which 

is  Gaussian  distributed  with  covariance  matrix  R,  =  I  + 

where  1  =  E(ninl)-  The  ergodic 
capacity  [4]  of  link  i  can  be  written  as 


Ci=E  log. 


1  —  /»,TT, .  Q.TTj  H  1 


where  the  expectation  is  taken  over  the  distribution  of  H. 

The  optimization  problem  to  maximize  the  system  capac¬ 
ity  can  be  defined  as 


L 

max  iCi  (la) 

i= 1 

subject  to  Tr(Qi)  <  1 
Q,>0. 


where  W;  are  pre-assigned  weights.  Due  to  the  nonconvexity 
of  this  problem,  only  a  suboptimal  solution  can  be  found. 


B.  Capacity  of  MIMO  Link  via  SVD 

Consider  the  MIMO  interference  system  with  a  rich  scat¬ 
tering  environment.  We  use  the  fact  that  a  MIMO  channel 
can  be  decomposed  into  a  number  of  parallel  independent 
channels.  By  multiplexing  the  independent  data  onto  these 
independent  channels,  we  can  increase  the  data  rate  signifi¬ 
cantly. 

Singular  value  decomposition  (SVD)  can  be  used  to  obtain 
the  independent  channels  for  a  MIMO  link.  Consider  a 
MIMO  link  with  nr  x  rit  channel  gain  matrix  H  that  is 
known  to  both  the  transmitter  and  the  receiver.  The  channel 
gain  matrix  H  is  decomposed  into 

H  =  UAVt, 

where  the  nr  x  nr  matrix  U  and  the  rit  x  rit  matrix  V  are 
unitary  matrices,  and  the  nr  x  nt  matrix  A  is  a  diagonal 
matrix  of  singular  values  {oy}  of  H.  Due  to  considering  a 
rich  scattering  environment,  the  matrix  H  is  full  rank  and 
Rank(  H)  =  M  =  min  (nr,nt). 

The  parallel  decomposition  of  a  MIMO  link  is  imple¬ 
mented  by  applying  transmit  precoding  at  the  transmitter  and 
receiver  shaping  at  the  receiver.  It  is  worth  pointing  out  that 
these  parallel  channels  are  independent  and  do  not  interfere 
with  each  other,  but  subject  to  the  total  power  constraint. 

The  channel  gain  matrix  H(J  from  the  transmitters  of 
eigenchannels  of  link  j  to  the  receivers  of  eigenchannels 
of  link  i  [i  f  j)  shown  in  Fig  1  can  be  written  as 

Hij  =  UjHijVj,  (2) 

where  U,  and  V,  are  from  the  singular  value  decomposition 
of  the  channel  gain  matrices  Hj  j  and  Hj  ;.  II,  defined  in 
(2)  is  verified  in  the  proof  of  Theorem  1.  Then,  Hlfk-  denotes 
the  channel  gain  from  the  transmitter  of  the  kth  eigenchannel 
of  link  j  to  the  receiver  of  the  Ith  eigenchannel  of  link  i. 

Theorem  1:  By  applying  SVD  decomposition  for  L  links 
in  the  MIMO  system  with  interference,  the  ergodic  capacity 
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Fig.  1.  The  channel  gain  matrix  from  the  transmitters  of  eigenchannels 
of  link  j  to  the  receivers  of  eigenchannels  of  link  i 


of  link  i  can  be  written  as 


Q=  El  log2 


PiAiQjxAj 


I  +  E.|  ftj-Hij QjxH^. 


where  A,  is  a  diagonal  matrix  of  singular  values  of  Hj.j,  Q,* 
is  the  covariance  matrix  of  transit  signal  x  before  precoding 
which  is  a  diagonal  matrix  diag{P.j,  ...  PiM\  where  Pik  is 
the  transmit  power  at  the  kth  eigenchannel  of  link  i,  and  the 
expectation  is  taken  over  the  distribution  of  H.  Furthermore, 
since  the  matrix  Q,*  of  link  j  is  a  diagonal  matrix  diag[Pj1 


...  P 


3M  J 


Hence, 


E 


Ci  = 

( 

log2 

V 


{)  A  Q.x  A 


t 


I  +  E4,  Pi, 3  Efcli  Pj>  (hm)  fc  (hm) 


where  ( H 


is  column  k  of  H 


The  proof  can  be  found  in  appendix.  It  is  easy  to  verify 

that  Ri  =  I  +  EjLi Pij  Efeli  Pjk  (h,^  *  is 

a  Hermitian  matrix  by  checking  R,  =  r].  By  applying  the 
fact  that  the  Hermitian  matrix  AA  is  positive  semi-definite 
for  any  matrix  A,  we  can  show  that  R,  is  a  positive  definite 
Hermitian  matrix. 

The  entries  of  H ij  are  i.i.d  complex  Gaussian  random 
variables  with  zero  mean  and  unit  variance  because  the 
unitary  matrices  Uj  and  Vj  do  not  change  the  distribution 
of  entries  of  H,  r  Hence,  the  non-diagonal  elements  of 


H, 


\  t ' 

H ij  J  have  zero  mean  and  the  diagonal  ele¬ 
ments  have  unit  mean.  For  large  L  and  M,  the  matrix 

C  =  Ej=i,j/i  Pi ,j  Efc=i  Pjk  k  can 


ap¬ 
proximated  by  a  diagonal  matrix  with  the  diagonal  elements 


of  C. 

The  aggregate  interference-plus-noise  at  the  receiver  of 
the  Ith  eigenchannel  of  link  i  from  the  eigenchannels  of 
other  links  has  the  covariance  C,  which  can  be  written  as 

L  M  , 

c*,«i+  E  ■ 

k=  1 

Then,  the  approximately  ergodic  capacity  of  link  i  can  be 
written  as 


Ci 


E 


PivlPj^ 


(3) 


where  the  expectation  is  taken  over  the  distribution  of  H, 
is  the  singular  value  of  the  Ith  eigenchannel  of  link  i,  and 
Pit  is  the  transmit  power  at  the  Ith  eigenchannel  of  link  i. 

In  summary,  the  channel  covariance  coefficient  DXrn  Vji 
from  the  transmitter  of  the  nth  eigenchannel  of  link  y  to  the 
receiver  of  the  mth  eigenchannel  of  link  x  can  be  represented 


as 

(  ctL 

if  x  =  y  and  m  =  n 

Dxm,yn  =  \  0  f 

if  x  =  y  and  m  f  n 

1  T/m , n  f  jErm,n  \ 

if  x  7 -  y. 

C.  The  Set  of  Discrete  Data  Rates  and  Powers 

Suppose  that  there  are  K  modulation/coding  schemes  and 
S  transmit  powers  available  for  the  mth  eigenchannel  of 
link  x ,  then  associated  with  each  MIMO  link  x  is  the  set 
of  logical  links  xm^k,s  with  1  <  m  <  M,  1  <  k  <  K  and 
1  <  s  <  S.  We  define  an  assignment  v  as  a  set  of  logical 
links  that  can  transmit  simultaneously.  In  other  words,  the 
assignment  specifies  which  links  are  transmitting,  which 
eigenchannels  are  active,  their  bit-rates,  and  their  transmit 
powers.  Specifically,  an  assignment  v  6  {0,  l}ix Mx A xS, 
where  vXm  k  s  =1  implies  that  the  MIMO  link  x  is 
transmitting  at  the  mth  eigenchannel,  the  kth  bit-rate  and 
the  sth  normalized  power. 

The  set  of  data  rates  are  selected  from  the  noninterference 
case.  The  transmitters  of  each  user  always  assume  that  there 
is  no  interference  from  other  users.  For  a  particular  SNR  = 
px ,  the  mth  eigenchannel  of  link  x  can  support  a  maximum 
data  rate  up  to 


CXm  —  l0g2  (l  +  PxaxmPmax)  ; 


where  Pmax  =  1  is  the  normalized  maximum  power. 
However,  there  exists  interference  when  there  are  more 
than  one  link  transmitting  simultaneously,  and  the  data  rate 
CXm  cannot  be  achieved. 

In  order  to  accommodate  the  interference  from  other  users, 
let  us  define  /?  as  the  buffer  size  to  reduce  the  sensitivity  to 
interference.  Therefore,  the  maximum  data  rate  of  the  mth 
eigenchannel  of  link  x  is 

CZX  =  1°§2  (i  + 

where  /3  =  2 dB  is  used  in  this  work. 
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Once  the  maximum  data  rate  of  eigenchannel  link  xm  is 
obtained,  a  wide  range  of  possible  approaches  are  available 
for  selecting  the  set  of  data  rates.  The  obvious  option  is  to 
select  K  linearly  separated  data  rates  up  to  (7™ax  denoted  by 
DR=[C™* / K  ...  C™^].  Then,  the  corresponding  SINK 
threshold  to  achieve  the  kth  data  rate  can  be  represented  as 

SINRf  =  2kxC*™/K  -  1. 

'Lm,k 

In  order  to  achieve  the  desired  data  rate,  the  SINR  at 
the  receiver  of  the  eigenchannel  must  be  no  less  than  the 
corresponding  SINR  threshold. 

Similarly,  the  set  of  powers  can  be  selected  as  S  linearly 
separated  transmit  powers  denoted  by  [1/5  ...  1].  The  num¬ 
ber  of  logical  link  for  each  MIMO  link  is  up  to  M  x  K  x  S, 
and  the  total  number  of  logical  links  in  the  ad  hoc  network 
is  up  to  L  x  M  x  K  x  S  which  detennines  the  computation 
complexity  of  the  optimization  problem  defined  in  (5). 

D.  Integer  Programming  Based  Algorithm 

Let  RXm  k  s  be  the  data  rate  over  the  logical  link  xm^,s, 
and  SINR (  denotes  the  SINR  threshold  required  to 
achieve  the  kth  data  rate  over  the  mth  eigenchannel  of 
link  x.  Note  that  SINR k  depends  on  the  link  x.  The 
aggregate  interference-plus-noise  at  the  receiver  of  the  mth 
eigenchannel  of  link  x  has  the  covariance  Qx  which  can  be 
written  as 

L  M  K  S 

Cxm  =  X^  Px.y  X!  Dxm,yn  X!  X!  Pyn,t,Qvyn.t.Q  +  1) 

y—l,y^x  n—  1  t=  1  q= 1 
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Fig.  2.  Capacity  per  User  for  nt  =  nr  =  2  and  L  =  2. 

where 

=SINRihmkx 

L  M  K  S  \ 

X  Px’V  X  Dxrn.,Vn  X  X  PVn,t,q  +  1  I  ' 

y—l^y^x  n=  1  t=  1  q— 1  J 

Constraint  (5c)  shows  that  the  total  power  at  the  transmitters 
of  link  x  must  be  no  more  than  the  normalized  power  1. 
Constraint  (5d)  enforces  that  only  one  logical  link  can  be 
active  for  each  eigenchannel  of  link  x. 

IV.  Numerical  Experiments 


where  PVri  t  is  the  transmit  power  of  logical  link  yn,t,q- 
The  logical  link  xm^,s  t°  be  active  must  satisfy 

PxDx  xmPXmik'S  >  SINRth^ 

Xm 

which  is  equivalent  to 


PxDxm,xmPxm,k,s  -  SINRfm  kCXm  >  oo  x  (vXm  k  e  -  1)  . 

(4) 

Since,  when  vXm  k  =  0,  (4)  is  always  true  (assuming  that 
oo  x  0  =  0).  Of  course,  in  practice,  oo  is  replaced  with 
some  large  number.  While  further  tuning  is  possible,  we  have 
found  that  Problem  (1)  can  be  approximated  by  solving  the 
following  integer  programming  problem. 


L  M  K  S 

“f  X  X  X  X  (5a) 

x—1  m—  1  k= 1  s=l 

subject  to 

PxDxm,xmPxm,k,a  -  SINR?m  k CXm  >  I-,,,.  ,.,.  - 

for  each  logical  link  xm^jS  (5b) 


M  K  S 

X  X  X  Px™,k,svxm,k,B  <  1  for  each  link  x  (5c) 

m=  1  k=  1  s=  1 
K  S 

X^  X/ Vxm  k  -  —  ^  ^or  eac^*  e*gencbannel  of  link  x  (5d) 

k= 1  s=l 


Vxm,k:e  G  {0,  1}  I 


In  order  to  compare  the  performance  of  different  ap¬ 
proaches,  i.e.,  GP,  QN  and  the  proposed  IP  method,  a 
large  number  of  simulations  are  performed.  The  results  are 
averaged  on  high  number  of  randomly  generated  channel 
matrices.  For  simplicity,  we  always  use  nt  nr  -  2,  and 
the  number  of  data  rates  and  powers  for  each  eigenchannel 
are  K  =  8  and  5  =  4,  respectively.  CPLEX  vlO  is  used  to 
solve  the  integer  programming  problem. 

A.  Capacity  Per  User 

In  the  symmetric  case,  SNR  and  INR  are  the  same  for 
all  users  where  =  p  and  pi  ;j  =  p  for  all  i  and  j. 

As  shown  in  [1],  [2],  the  ergodic  mutual  information  is 
relatively  independent  to  the  different  initial  conditions.  The 
GP  and  QN  methods  converge  in  most  cases,  and  the  scare 
nonconvergence  cases  are  excluded  in  the  simulations. 

Figs  (2-3)  show  the  ergodic  mutual  information  per  user 
versus  INR)?/)  for  different  SNR  values  for  two  and  five 
users,  respectively.  Capacity  per  user  decreases  with  the 
increase  of  INR,  the  decrease  of  SNR,  and  the  increase  of 
the  number  of  users. 

Not  surprisingly,  GP  and  QN  methods  achieve  better 
performance  than  IP  method  when  INR  is  small.  [1]  has 
shown  that  the  ergodic  mutual  information  in  the  network 
is  a  convex  function  of  links’  covariance  matrices  when 
INR  is  sufficiently  small.  Hence,  both  GP  and  QN  methods 
can  achieve  the  optimal  performance.  However,  IP  method 
achieves  suboptimal  solution  due  to  several  approximations 
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Fig.  3.  Capacity  per  User  for  rat  =  nr  =  2  and  L  =  5. 


such  as  the  discrete  set  of  data  rates,  discrete  powers,  and 
the  interference  margin  /3.  As  shown  in  the  Figs,  the  results 
of  IP  method  are  still  very  close  to  the  optimal  solution. 

In  the  case  where  1NR  is  large,  IP  method  achieves  better 
performance  than  GP  and  QN  methods.  The  reason  is  that 
GP  and  QN  methods  solve  the  optimization  problem  with 
the  constraint  Tr(Q)  =  1.  In  other  word,  all  links  transmit 
with  full  power,  which  is  obviously  far  from  optimal  for  the 
strong  interference  case.  In  fact,  the  optimization  problem 
should  be  solved  with  the  constraint  Tr{Q)  <  1,  which 
means  scheduling  is  necessary  and  some  links  are  selected 
to  transmit  when  the  INR  is  large. 

B.  Average  Number  of  Active  Links 

Figs  4  shows  the  average  number  of  active  links  obtained 
by  IP  method  for  five  users.  Fixing  the  value  of  SNR,  the 
average  number  of  active  links  decreases  with  the  increase 
of  INR.  On  the  other  hand,  fixing  the  value  of  INR,  the 
average  number  of  active  links  increases  with  the  increase 
of  SNR  except  the  low  INR  and  high  INR  cases  where  all 
links  transmit  or  only  one  link  is  selected  to  transmit. 

When  INR  is  sufficiently  small,  all  users  transmit  simul¬ 
taneously,  which  can  be  verified  by  the  case  INR=—  5dB. 
When  INR  is  large  such  as  INR=20dB,  the  number  of  active 
users  is  close  to  one.  Therefore,  transmitting  one  link  at  a 


L  =  5. 

time  is  better  than  transmitting  all  links  simultaneously  with 
full  power  for  high  INR  case.  Generally  speaking,  it  is  better 
to  select  links  to  transmit  when  INR  is  large. 

C.  Average  Power  of  Active  Links 

Figs  5  shows  the  average  power  per  active  link  obtained  by 
IP  method  for  five  users.  As  we  know,  GP  and  QN  methods 
assume  that  all  links  are  active  and  transmit  with  frill  power. 
As  shown  in  Figs  (4-5)  for  IP  method,  all  links  transmit 
with  full  power  for  the  weak  interference  case,  and  only 
one  link  is  selected  to  transmit  with  full  power  for  the  strong 
interference. 

The  scenarios  between  these  two  extreme  cases  are  com¬ 
plicated.  As  we  can  see  from  Fig  5,  the  average  power  per 
active  link  is  lower  for  larger  SNR  case.  The  reason  is  that 
more  links  are  selected  to  transmit,  and  it  is  necessary  to 
decrease  the  transmit  power  in  some  links  to  reduce  the 
interference  to  other  active  links,  resulting  in  lower  average 
power  per  active  link. 

D.  Average  Power  of  Active  Eigenchannel  Links 

Fig  6  shows  the  average  power  per  active  eigenchannel 
obtained  by  IP  method  for  five  users.  The  transmit  power 
per  eigenchannel  is  close  to  0.5  for  weak  interference  case 
{SNR  =  20 dB  and  INR  =  —5 dB),  which  matches  the 
result  in  [9].  Whereas,  the  transmission  power  per  eigen¬ 
channel  is  also  close  to  0.5  for  strong  interference  case 


5 


(INR  =  20 dB),  which  is  not  in  conflict  with  the  result  in 
[9].  The  reason  is  that  only  single  link  is  selected  to  transmit 
for  strong  interference,  which  indicates  that  the  active  link 
transmits  without  interference. 

V.  Conclusion 

In  this  paper,  we  study  the  capacity  optimization  problem 
of  MIMO  links  with  interference.  Due  to  the  nonconvexity 
of  the  capacity  problem,  some  suboptimal  solutions  have 
been  proposed  in  the  previous  works,  such  as  a  Gradient 
Projection  method  [1]  and  a  Quasi-Newton  method  [2].  Both 
of  these  methods  assume  Tr{Q )  =  1,  implying  always 
transmitting  with  full  power. 

We  derive  the  capacity  for  MIMO  links  suffered  from 
cochannel  interference  where  each  MIMO  link  is  decom¬ 
posed  via  SVD,  and  present  an  Integer  Programming  based 
algorithm  to  solve  the  capacity  optimization  problem.  Com¬ 
pared  with  GP  and  QN  methods,  IP  method  achieves  lower 
capacity  for  the  case  of  weak  interference  and  higher  ca¬ 
pacity  for  the  case  of  strong  interference.  In  the  case  of 
weak  interference,  GP  and  QN  methods  achieve  the  optimal 
solution  due  to  the  convexity  of  the  optimization  problem. 
In  the  case  of  strong  interference,  transmitting  one  link  at  a 
time  is  better  than  transmitting  all  links  simultaneously  with 
full  power.  Generally  speaking,  scheduling  links  to  transmit 
is  more  efficient  for  the  case  of  strong  interference. 

Appendix 


The  channel  gain  matrix  Hj j  defined  in  (2)  is  verified  and 
Hi, j  =  UtHjjVj.  Then, 


7  = 


f>,  A ,  Q  A 


We  know  that  the  matrix  Qj*  of  link  j  is  a  diagonal  matrix 
diag[Pj1  ...  PjM]  where  Pjk  is  the  transmit  power  at  the  kth 
eigenchannel  of  link  j.  Hence, 


7  = 


1  + 


/ 1  7  (  j .  x  7 


I  +  Ej=l,j^i  Pi,j  Efc= 1  Pjk  (Hi,j^  (HijJ 


where  [HjjJ  is  column  k  of  H,j.  Therefore,  the  ergodic 
capacity  of  lini:  i  can  be  written  as 


Ci  =  E  (log2  7) 


where  the  expectation  is  taken  over  the  distribution  of  H.  ■ 
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Proof  of  Theorem  1:  Applying  the  transmit  precoding 
at  link  i,  the  covariance  matrix  of  the  transmit  data  x  can 
be  written  as 


Q  =  V  Q  xV 


t 

i  ’ 


where  x  =  Vx.  Then, 
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i  +  p  ii  q  n;  r,  1 

I  +  Pi(UiAlvJ)(V,Qiiv])(VlAjuJ)R- 
I  +  ftUiAiQisAtutRr1!. 


1 


Applying  the  fact  that 

|Im  +  AB|  =  |I„  +  BA| 


where  A  is  a  m  x  n  matrix  and  B  is  a  n  x  m  matrix. 
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I  +  PjAiQj*  A-U-R^Uj 
I  +  PiAlQiiAt(ujRlUi)-1 


Because  of  the  precoding  of  link  j,  we  get 
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